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ABSTRACT 

Context. RU Lup is a Classical T Tauri star with unusually strong emission lines, which has been interpreted as manifestations of accretion. 
Recently, evidence has accumulated that this star might have a variable radial velocity. 

Aims. We intended to investigate in more detail the possible variability in radial velocity using a set of 68 high-resolution spectra taken at the 
VLT (UVES), the AAT (UCLES) and the CTIO (echelle). 

Methods. Using standard cross-correlation techniques, we determined the radial velocity of RU Lup. We analysed these results with Phase- 
dispersion minimization and the Lomb-Scargle periodogram and searched for possible periodicities in the obtained radial velocities. We also 
analysed changes in the absorption line shapes and the photometric variability of RU Lup. 

Results. Our analysis indicated that RU Lup exhibits variations in radial velocity with a periodicity of 3.71 days and an amplitude of 2.17 
km s _1 . These variations can be explained by the presence of large spots, or groups of spots, on the surface of RU Lup. We also considered a 
low-mass companion and stellar pulsations as alternative sources for these variations but found these to be unlikely. 
Conclusions. 
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1. Introduction 

T Tauri stars (TTS) show hydrogen emission lines superim- 
posed on relatively normal late-type photospheric spectra in the 
optical spectral region. A fraction of the TTS have very strong 
and broad emission lines, and in addition emission lines of neu- 
tral and once ionised metals, and in some cases He. These stars 
are nowadays referred to as the classical T Tauri stars (CTTS) 
since they were the first to be studied spectroscopically by 
Joy ( 1945). CTTS are identified with pre-main-sequence low- 
mass stars, and many of their observed characteristics can be 
explained within the framework of magnetospheric accretion 
(see for example, Hartmann 1998 ). According to this model, 
gas is ionised at the interior of a circumstellar, dusty disk, af- 
ter which the plasma is accreted along stellar dipole magnetic 
field lines, and eventually dumped under virtually free fall at 
high latitudes, creating a shocked region that heats the photo- 
sphere locally. The energy released in this hot accretion shock 
then manifests itself as excess continuous emission. The model 
of magnetospheric accretion has gained support through the de- 
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tection of relatively strong (several kG) magnetic fields around 
CTTS (Guenther et al. [19991 Johns-KruU et al. |1999al Valenti 
& Johns-Krull |20Q4l ). 

Among the CTTS, there is a small group of stars with ex- 
tremely rich emission line spectra and veiling, which on oc- 
casions can obliterate the photospheric spectrum completely. 
Some of these "enhanced" CTTS vary dramatically and ran- 
domly in brightness, line profiles and veiling on time scales 
as short as hours, which is often interpreted as being related 
to changes in the instantaneous accretion rate. Still, what the 
driving force is behind these strongly enhanced accretion sig- 
natures is unclear. 

Several of the enhanced CTTS have turned out to be close, 
double-lined spectroscopic binaries. In DQ Tau (Basri et al. 
[T9971 UZ Tau E (Martin et al.|2005]) and V4046 Sgr (Stempels 
& Gahm 2004), accretion signatures appear to be modulated 
by the binary orbit. Other "enhanced" CTTS show velocity 
changes of small amplitude in single photospheric lines and 
may be single-lined spectroscopic binaries (Melo 2003 ; Gahm 
2006 ). It is also possible that the changes in radial velocity are 
caused by starspots moving across the stellar disk. It may be 
difficult to discriminate between the two scenarios. An example 
is the extreme CTTS RW Aur A for which Petrov et al. (USB 
considered a low-mass binary companion. However, they also 
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found reasonable fits to the line variations by modeling spot- 
induced absorption line variations. 

Within this framework we will discuss the well-known and 
extremely active CTTS RU Lup, for which we recently discov- 
ered a possible periodicity in radial velocity (Gahm et al. 2005 ). 
RU Lup is well-studied both spectroscopically and photometri- 
cally, displaying very strong and variable accretion signatures. 
The properties of RU Lup were described in depth by Herczeg 
et al. d2505l and Stempels & Piskunov ( 120021 . 

2. Observations 

We observed RU Lup with the UVES spectrograph at the 8-m 
VLT/UT2 in Chile during three observing runs, April 15- 
16, 2000, April 15-18, 2002 and August 14-17, 2005. We 
used UVES in dichroic-mode, to obtain simultaneous expo- 
sures with the blue and red arm of the spectrograph. In this 
way, the wavelength coverage is almost continuous, and ranges 
from 3500 to 6700 A. High-quality spectra, normally at a res- 
olution of R « 60 000, could be obtained in exposures as short 
as 15-30 minutes with the blue arm (S/N-ratio ~ 50), and 10 
minutes with the red arm (S/N-ratio ~ 150). In total we ob- 
tained 29 blue and 66 red spectra. 

In addition to the UVES data, we also obtained one spec- 
trum of RU Lup taken with the echelle spectrograph at the 
4-m CTIO, on July 17, 1998, and one spectrum taken with the 
UCLES spectrograph at the 4-m AAT, on July 11, 2005. 

All data were reduced with the IDL-based package 
REDUCE (Piskunov & Valenti 2002). This package uses ad- 
vanced techniques for accurate bias-subtraction, modelling of 
scattered light and flat-fielding, and performs optimal order ex- 
traction through 2D modelling of the individual orders. 

3. Analysis 

3.1. A synthetic template spectrum for RU Lup 

The spectrum of RU Lup contains a large number of strong 
emission lines of e.g. H, Ca n H & K, Fe i and n, [O i] and 
[S n]. Although line emission and veiling contaminate large 
parts of the spectrum, there are shorter regions where narrow, 
photospheric absorption lines can be seen. In these regions we 
determined the photospheric contribution in RU Lup by cal- 
culating synthetic template spectra (see Stempels & Piskunov 
2002 2003). This method assumes that the "stellar" spectrum 
can be described by the sum of a featureless excess continuum 
and a normal late-type photospheric spectrum. The model at- 
mosphere, stellar parameters and atomic line data used are the 
same as in Stempels & Piskunov (2003). Other advantages of a 
synthetic template spectrum is that it allows to accurately deter- 
mine the stellar radial velocity and absorption line shapes from 
cross-correlation. 

3.2. Veiling 

RU Lup exhibits both photometric and spectroscopic varia- 
tions. An important contributor to this variability is the veiling 
continuum. Veiling manifests itself as a featureless continuum 



Table 1. Modified Julian Date, measured radial velocities and 
the Bisector Inverse Slope (BIS, see Sect. |3.4| ) for RU Lup. The 
error (cr) of our radial velocity measurements is 0.16 km s _1 . 
Values marked with * are taken from Melo (120031) . but trans- 
ferred to MJD. The observation marked with ** was obtained at 
the CTIO, and the observations marked with *** at the AAT. 
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that decreases the depth of photospheric absorption lines. The 
veiling factor is defined as the ratio of the pseudo-continuum 
and the stellar continuum. 

Part of the current dataset (the data from 2000) were used 
by Stempels & Piskunov (2002 ) for a detailed analysis of the 
veiling in RU Lup. We applied the same method on the en- 
tire dataset and found an average veiling factor of 3.4, but with 
strong short-term variability. Some observations of RU Lup ob- 
tained in 2005 have veiling factors as high as 9. 
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3.3. Radial velocities 

We used traditional cross-correlation with a (synthetic) tem- 
plate spectrum over the spectral region 6000 to 6030 A, which 
is entirely free from line emission. Here, the S/N-ratio of our 
observed spectra is ~ 125, and the ^-statistics of our measure- 
ments yield typical errors of cr — V27VV 8 0.16 km s- 1 . 
Correction for possible instrumental drifts was done by cross- 
correlating the telluric oxygen lines around 6300 A with the 
corresponding lines in the NSO solar spectrum (Kurucz et al. 
1984). The corrections turned out to be about +0.2 km s _1 
within a night. These corrections, as well as the heliocentric 
velocity correction, were applied to each observed spectrum. 
Finally, we obtained radial velocities for RU Lup by cross- 
correlating the (corrected) observed spectra with the synthetic 
template spectrum. 

The obtained velocities, given in TableQ] range from about 
-4.0 to +1.0 km s _1 , which clearly shows that RU Lup displays 
changes in radial velocity that are much larger than the velocity 
corrections due to instrumental drifts or any other uncertainty. 

We also checked the quality of our reduction method by 
performing a similar analysis on spectra of a radial velocity 
standard star. The scatter in radial velocity was less than 0.05 
km s _1 , indicating that the internal precision of our radial ve- 
locity measurements is satisfactory. 

3.4. Line bisectors 

The observed amplitude of v ra( j is ~ 2.2 km s . This is well 
within RU Lup's v sin i value of 9.0 km s _1 . It is thus important 
to consider whether these are genuine line shifts or possibly 
due to asymmetries within the line profile. 

At medium resolution or low v sin i it may be difficult to 
determine asymmetries from individual line profiles. Instead 
of using single lines, an "average" line profile can be obtained 
from the cross-correlation function (CCF) of the photospheric 
absorption lines. The bisector of the CCF then provides a mea- 
sure for the "average" line asymmetry. The method of calcu- 
lating CCF bisectors and its applications are described in more 
detail by Dall et al. ( 120061) . In our analysis we used the same 
CCF as the one we recovered from the radial velocity measure- 
ments in the previous section, except that we used a template 
without rotational broadening. Cross-correlation with such a 
narrow template enhances any asymmetries in the CCF. 

There exist several definitions of measures to quantify the 
asymmetry of a CCF bisector. In our analysis we will use the 
so-called Bisector Inverse Slope (BIS, see Queloz et al. 2001) 
defined as BIS - v, - v/,, where v, is the mean bisector ve- 
locity between 10% and 40% of the line depth, and Vj is the 
mean bisector velocity between 55% and 90% of the line depth. 
We found that that there is a detectable line asymmetry in RU 
Lup, and that the BIS varies between approximately +1 and 
-1 km s _1 , see TableQ] The relation between the BIS and the 
observed radial velocities is discussed in Sect. 15.11 



Phase-dispersion minimization 




0.0 0.1 0.2 0.3 0.4 0.5 

Frequency (cycles/day) 



20 z Mi Lomb-ScarglB periodagram 

~ iii of datawindow - 

-0.2 -0.1 0.0 0.1 0.2 

Frequency (cycles/day) 

Fig. 1. Upper panel: Comparison of the results of phase disper- 
sion minimization and the Lomb-Scargle periodograms based 
on the observed data from Table Q] Some of the possible so- 
lutions, with periods of 3.63293, 3.71058, 4.45633, 5.58098, 
and 7.94913 days are highlighted. For the Lomb-Scargle pe- 
riodogram, three levels of the false-alarm probability are in- 
dicated. Lower panel: The Lomb-Scargle periodogram of the 
data window function, with the peak at zero frequency re- 
moved. For purpose of clarity, a different frequency domain 
was used for this panel. 

4. Periodicity 

4.1. A 3.7-day period in the radial velocity 

The relatively smooth changes in radial velocity within a night 
prompted us to search for periodic radial velocity changes. We 
also searched the literature for published radial velocity mea- 
surements of RU Lup, and we included data published by Melo 
(2003) in our analysis (see TableQ]). Other published velocities 
of RU Lup lack information on exact times of observation. 

The combined dataset we have at our disposal is inhomo- 
geneously sampled in time, sometimes with several observa- 
tions in short periods of time. The differences between the ra- 
dial velocities of consecutive spectra is small. Therefore, in or- 
der to focus on the long-term evolution of the radial velocity, 
rather than to contaminate the period analysis with the very 
short-term variations, we averaged all consecutively observed 
observations, leaving us with a set of 38 radial velocity mea- 
surements, including the 7 values published by Melo. For our 
analysis we assumed an uncertainty in the radial velocities of 
cr = +1 km s _1 . 

We analysed the data using two different period-detection 
methods. These methods are the Lomb-Scargle periodogram 
(according to the prescription of Home & Baliunas 1986 ) and 
phase-dispersion minimization (Stellingwer fl 19781) . The result- 
ing diagrams for these methods are shown in Fig.Q] 

An important feature of the Lomb-Scargle periodogram is 
that the height of the peaks in the periodogram allows to make 
an estimate of the false-alarm probability (FAP). We performed 
Monte-Carlo simulations to determine the probability that, for 
a given signal in the periodogram, normally-distributed random 
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Fig. 2. Upper panel: All observations of RU Lup, folded with 
the 3.71058 day period. Observations marked with + are taken 
at the VLT, with * at the AAT and A at the CTIO. Observations 
by Melo (2003 ) are marked with □. Lower panel: The Lomb- 
Scargle periodograrn after pre-whitening of the data with the 
3.71058 day period. Levels of the false-alarm probability are 
indicated. 

fluctuations give rise to a false detection of periodicity (i.e. the 
FAP). In Fig. Q] we indicate the levels in the periodograrn cor- 
responding to a FAP of 0.05, 0.01 and 0.001. The peaks in our 
periodograrn clearly exceed these levels, indicating that there 
is detectable periodicity in our data, although it is not directly 
clear which period fits the data best. 

The difficulty in determining the correct period with the 
Lomb-Scargle periodograrn is for a large part due to the inho- 
mogeneous distribution of our observations in time, producing 
a broad window function with strong sidelobes. (See Fig. [1). 
Also, the Lomb-Scargle periodograrn relies on the implicit as- 
sumption that the signal is purely sinusoidal. Non-sinusoidal 
signals might be detected, but with lower peaks of significance. 
Phase-dispersion minimization is unbiased towards the shape 
of the radial velocity curve, and might detect signals that are 
difficult to detect with the Lomb-Scargle periodograrn. 

As can be seen from Fig.Q] the Lomb-Scargle periodograrn 
and Phase-dispersion minimization gave a number of possible 
solutions. Presented with more than one clear solution, we per- 
formed pre-whitening of the radial velocity data and recalcu- 
lated the periodograms. We found that the pre-whitened peri- 
odograrn of all periods contained significant levels of power, 
except for the solution with P = 3.71058 days. The pre- 
whitened periodograrn for this period is shown in the lower 
panel of Fig. [2] The remaining power is clearly less than the 
lowest levels of the FAP. After pre-whitening, we can also es- 
timate the uncertainty in the period (Home & Baliunas, 1986) 
finding that P = 3.71058 + 0.0004 days. 

Although this period is well identified by Phase-dispersion 
minimization, it does not produce a strong peak in the 
Lomb-Scargle periodograrn. However, recalculating the Lomb- 
Scargle periodograrn using calculated radial velocities pro- 
duced a diagram very similar to the periodograrn of the ob- 
served data, including peaks of other candidate periods found 



in the original Lomb-Scargle periodograrn This supports our 
conclusion that the most likely period that fits our data has 
P = 3.71058 days. 

If the period derived above is identical to the rotational pe- 
riod, this allows us to determine the inclination of RU Lup. 
Herczeg et al. (2005) derived a radius of 1 .64 R0 and a mass 
of 0.65 Mq from the stellar luminosity, and if P ml * 3.71 days 
and v sin / = 9.0 + 0.9 km s _1 (Stempels & Piskunov 2002) then 
i » 24°. This is compatible with other observational arguments 
(Herczeg et al. 2005 , Stempels & Piskunov 2002 ) that indicate 
that RU Lup is observed close to pole-on. 

4.2. Photometric periodicity? 

RU Lup is well-known to exhibit strong erratic photometric 
variability. Early on, Hoffmeister ( 1965 ) and Plagemann (1969 ) 
reported a photometric period of 3.7 days. Giovannelli et al. 
(1991) found a weak signal corresponding to 3.7 days in the 
power spectrum when re-analysing the early data, and Drissen 
et al. (119891) found some evidence of a 3.7 day period in their 
linear polarization measurements. However, later analyses in- 
cluding newer data indicated that the 3.7 day photometric pe- 
riod is not significant (Gahm et al. [1993; Giovannelli et al. 
[19941 

From 2001 until 2005, RU Lup was also monitored in 
the V-band by the All Sky Automated Survey (ASAS, see 
Pojmanski 2002 ). This period of time covers most of our spec- 
troscopic observations of RU Lup. We have analysed the ASAS 
dataset of RU Lup, obtaining K = 11.23 ± 0.26, but no period- 
icity. Neither did we find any short-term correlation between 
the measured radial velocities in individual observing runs and 
simultaneous photometric data. We performed several numeri- 
cal tests on the ASAS data and found that the semi-amplitude 
of any possible photometric periodicity must be less than 0.12 
mag. Periodic variations larger than this would have been de- 
tected from these data with high probability. 

5. Discussion 

Periodic radial velocity changes in photospheric absorption 
lines may be generated in different ways, and below we will 
discuss possible sources such as rotational modulation of dark 
and/or bright spots on the stellar surface, binarity and pulsa- 
tions. 

5.1. Spots 

Spots can distort the shape of photospheric line profiles, an 
effect exploited in Doppler Imaging. Unfortunately, RU Lup 
has a vsin/ of only 9 km s _1 making it unsuitable as a target 
for Doppler Imaging. Still, spots may be an important candi- 
date to consider as a cause of the observed velocity changes in 
RU Lup. 

A critical test of the effects of spots on the stellar line pro- 
files can be performed using line bisector analysis. An inhomo- 
geneous distribution of spots on the stellar surface will make 
the line profiles appear asymmetric, which is detectable as a 
change in the line bisector. In addition to changes to the line 
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Fig. 3. Upper panel — left: The relation between the Bisector Inverse Slope (BIS, see text) and the observed radial velocity of 
RU Lup. The slope of the best least-squares fit is shown. Upper panel - right: The observed radial velocity (+) and the BIS (O) 
folded with the 3.71058 day period. The values of the BIS have been plotted with respect to the mean average radial velocity. 
Lower panel: As in the upper panel, but for a set of synthetic spectra of a star with a photosphere of T — 4000 K and a spot of 
T = 3400 K, a latitude of 60° and a radius of 35°. 



bisector, spots may cause variations in the apparent radial ve- 
locity of the star, within a fraction of its v sin i value. 

Any correlation between the BIS and the apparent radial ve- 
locity is often used as a strong indication that the shape of the 
line profiles changes are due to some process other than gen- 
uine shifts in radial velocity. For example, Queloz et al. ((2001) 
used an anti-correlation between the BIS and the observed ra- 
dial velocity to show that spots can cause a periodic signal in 
the apparent radial velocity. Similarly, Santos et al. (2002 ) used 
a positive correlation between the BIS and the observed radial 
velocity to argue that faint light from an undetected compan- 
ion changes the intrinsic line profile. Dall et al. (20© also 
found a positive correlation for EK Eri, and related this to sur- 
face activity, or possibly an unseen companion. The absence 
of a correlation is normally interpreted as a strong argument 
for genuine changes in radial velocity due to a companion (for 
example Santos et al. 120031 . 

In the top-left panel of Fig. [3] we show the relation between 
the BIS and the radial velocity we obtained for RU Lup. As 



one can see there is an anti-correlation between the BIS and 
the apparent radial velocity of RU Lup. The anti-correlation is 
more evident in the top-right panel, where we show both the 
BIS and the radial velocity as a function of phase. 

Although the anti-correlation between the BIS and the ra- 
dial velocity is a strong indication that surface structures might 
be responsible for the observed changes in radial velocity, it 
is important to quantify this effect in terms of spot size and 
changes in brightness. Therefore, we simulated the effect of 
starspots on the shape of photospheric lines using a simple 
model of a spotted star (see Hatzes 2002 ). 

In this model we assume parameters typical for RU Lup, 
r eff = 4000 K, vsinz' = 9 km s" 1 and i = 24°. The photo- 
sphere was represented by a synthetic K7 V spectrum, while for 
the spots we used synthetic spectra across a range of tempera- 
tures (3000 - 3800 K). The photospheric and spot spectra were 
then combined to construct disk-integrated synthetic spectra 
of a spotted star for several rotational phases and spot radii 
(0 - 40°). The resulting spectra cover the spectral region 6000 
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Fig. 4. Calculated effect of a stellar spot located at a latitude 
of 60° on the amplitude of radial velocity variations for a range 
of temperatures and spot radii. The dotted line indicates the 
observed amplitude of variations in RU Lup. 

- 6050 A and have a spectral resolution of about 60 000. For 
each spectrum we determined the corresponding radial veloc- 
ity v'rad, as well as the brightness and the bisector shape through 
cross-correlation with the spectrum of an unspotted K7 V star. 

From our model we find that one single spot at a latitude 
of =s 60° and on an otherwise undisturbed stellar surface will 
cause the most dominating effect (i.e. on radial velocity, bright- 
ness and line asymmetry). Although using one single spot is 
the simplest possible model to reflect the observed variations, 
this does not necessarily reflect the actual distribution. Doppler 
Imaging of other CTTS has shown that a distribution with sev- 
eral (groups of) spots is more realistic. We therefore interpret 
the results of our simulation only as a first-order indication of 
the effect of an asymmetric distribution of spots on the ob- 
served line profile. 

The shift in observed radial velocity due to a spot is a func- 
tion of both the spot size and the temperature contrast between 
the spot and the stellar photosphere, as is illustrated in Fig. [4] 
This figure shows that, in order to explain the observed am- 
plitude in radial velocity, the surface of RU Lup must have a 
relatively large spot coverage (similar to a single spot with a 
radius of > 30°), and that these spots can have (average) sur- 
face temperatures in the range 3000 - 3600 K. 

We did similar calculations for hot spots, assuming the 
spot to emit a featureless spectrum (black-body continuum) of 
7000 K (Calvet & Gullbring[l998]l. Because hot spots have a 
much larger surface brightness than the rest of the photosphere, 
their contribution to the total spectrum is much stronger than is 
the case for cool spots. We found that, although asymmetric 
hot spot may explain changes in radial velocity, the implied 
spot radius (^ 40°) would cause very strong levels of veiling 
(^ 10), which would correlate strongly with phase. Hot spots 
would also cause strong phase-related changes in the blue part 
of the spectrum. We found no evidence for either effect in our 
observations of RU Lup. In addition, CTTS tend to have typ- 
ical values of hot spot filling factors around 0.1 - 5% (Calvet 



& Gullbring[T998l Calvet et al. l2004l and not 25 - 50%. Non- 
polar hot spots also are transient phenomena related to (vari- 
able) accretion. We find it therefore unlikely that these short- 
lived, hot surface features could be responsible for the long- 
term velocity variations in RU Lup. This, of course, does not 
rule out active accretion and small hot spots on RU Lup. 

Our simulations also allowed us to estimate the magnitude 
of brightness variations due to spots. Although a cool spot with 
a radius of 40° will induce a relatively strong modulation in 
brightness, it will be difficult to detect this effect from photom- 
etry alone, since high levels of veiling will strongly reduce the 
brightness contrast of such variations. This reduction in bright- 
ness contrast between spot and star is also the reason why we 
were unsuccessful in detecting other spectroscopic signatures 
of spots, such as molecular bands from TiO. 

For each set of simulated spectra we determined the behav- 
ior of the BIS, using the same method and template as for the 
observed spectra. We found that the BIS of the simulated spec- 
tra shows an anti-correlation with rotational phase, and that the 
slope of this correlation depends weakly on both the spot size 
and temperature. From our set of models, a star with a 35° spot 
of T — 3400 K shows the best agreement with the observations, 
fitting the slope of the anti-correlation, as well as the shape and 
amplitude of the velocity curve. This is illustrated in the lower 
panels of Fig. [3] 

We conclude that it is possible to explain the observed 
changes in radial velocity in RU Lup with cool spots, and that 
this requires relatively large surface areas covered with spots. 
The consistency of the radial velocity variations, which allowed 
us to establish the rotational period, indicates a long lifetime 
for the spot structures on RU Lup. Other CTTS for which 
approximate spot sizes have been determined seem to follow 
this pattern. Doppler Imaging studies of Sz 68 (Johns-Krull & 
Hatzes [T997j l, DF Tau (Unruh et al. fT998l >, SU Aur (Petrov et al. 
[T9961 Unruh et al.|2004]and MN Lup (Strassmeier et al. l2003T l 
all recover large surface areas covered by spots. Modelling of 
the spot coverage fraction in active stars using molecular band 
analysis suggests that Doppler Imaging, and thus also bisector 
analysis, may underestimate the true spot filling factor due to 
features that do not change the line profile (O'Neal et al. 2004 ). 

The large spot filling factor we find for RU Lup could also 
be in agreement with the empirical relation between spot filling 
factor, photospheric temperatures and magnetic field strengths 
on active dwarfs derived by Berdyugina (2006). This relation 
would imply a strong magnetic field of the order of 3 kG in this 
star. 

5.2. Binarity 

One alternative explanation of the radial velocity variations in 
RU Lup may be that a low-mass secondary component orbits 
the star. An elliptic binary orbit with P = 3.71058 days and 
e - 0.41 is compatible with the periodic radial velocity varia- 
tions observed in RU Lup (see Fig. [2j. Both the amplitude and 
the shape of the velocity curve match the effect of a low-mass 
secondary orbiting RU Lup. 
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The mass of the secondary can be determined under the 
assumption that the axis of rotation are lined up, and that 
the system rotates synchronously. This is a normal assump- 
tion for close binary systems, because the time-scale for syn- 
chronization for such systems is relatively small (see Hilditch 
2001 and references therein). The mass function for a single- 
lined spectroscopic binary is f(m) - sin 3 i/(rti2 + mi) 2 = 
1.0361 • 10~ 7 (1 - e 2 f> 2 K\P, where e is the ellipticity, K x is the 
amplitude of the velocity variations of the primary in km s , 
and P the period of the system in days. Using m.\ = 0.65 M©, 
Ki = 2.17 km s _1 and i = 24°, we derive a secondary 
mass of 1112 = 0.027 M©, which is well within the brown 
dwarf regime. The semi-major axis of the primary follows from 
fli sin; = 1.9758 ■ 10~ 2 (1 - e 2 ) 1/2 ffiP giving a x = 0.36 R and 
«2 = m\lni2 ■ a\ — 8.5 Rq. 

Although the orbital solution fits the observed radial veloc- 
ity of RU Lup, we find the solution of a brown dwarf secondary 
unlikely. The obtained solution does not allow us to reproduce 
the observed correlation between the BIS and the radial veloc- 
ity. While a brown-dwarf binary seems unlikely, our observa- 
tions cannot dismiss a possible RS CVn-like scenario in which 
spots and a binary companion cause simultaneous and phased 
modulation of the radial velocity. 

5.3. Pulsations 

The loci of fundamental pulsation periods of contracting stars 
over the HR-diagram were calculated by Gahm & Liseau 
(1988). The position of RU Lup in the HR-diagram indi- 
cates possible periods on time-scales of hours rather than days 
(Gahm 2006). Unless there is a not yet identified mode of sur- 
face oscillations, it appears that pulsation is not the agent of 
periodic radial velocity changes. 

6. Conclusions 

Using standard cross-correlation techniques we have found ra- 
dial velocity variations in very high resolution spectra of the 
CTTS RU Lup. The variations span a small range of velocities 
(-4 to +0.5 km s _1 ) and appear to be periodic (P * 3.71 days). 

These radial velocity variations correlate with the slope of 
the bisector, indicating that they are probably related to spots 
on the stellar surface. We quantified this effect by simulating 
these variations with a model of a spotted star. We find that a 
large part of the surface of RU Lup has to be covered with long- 
lived spots of T ~ 3400 K to explain the observed changes in 
the bisector and radial velocity. A large spot coverage indicates 
that a surface magnetic field of about 3kG can be expected in 
RU Lup. 

We also considered binarity and pulsations as other possi- 
ble sources of radial velocity variations, but these turned out to 
be unlikely candidates. 

The 3.71 day period in radial velocity we recovered is sim- 
ilar to what was reported in some very early photometric stud- 
ies of RU Lup. Even though later investigations, including this 
work, do not confirm any photometric periodicity, it would be 
of value to repeat multi-colour monitoring, especially in the 



blue part of the spectrum, to search for underlying weak sig- 
nals of periodic nature. 
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